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Numerical calculations based on the recursive Green's functions method in the tight-binding ap- 
proximation are performed to calculate the dimensionless conductance g in disordered graphene 
nanoribbons with Gaussian scatterers. The influence of the transition from short- to long-ranged 
disorder on g is studied as well as its effects on the formation of a perfectly conducting channel. 
We also investigate the dependence of electronic energy on the perfectly conducting channel. We 
propose and calculate a backscattering estimative in order to establish the connection between the 
perfectly conducting channel (with g = 1) and the amount of intervalley scattering. 



I. INTRODUCTION 

The remarkable electronic transport properties of 
graphene have motivated numerous experimental and 
theoretical studiesji"— Of particular interest is the pos- 
sibility of fabricating narrow width graphene samples, 
called graphene nanoribbons (GNRs). By engineering 
the lateral confinement one can, in principle, create an 
electronic energy gap leading to a semiconductor behav- 
ior that allows for the development of novel electronic 
nanodevices and applications/^ 

The observation of conductance quantization in GNRs 
turned out to be more difficult than anticipated^. The 
reason is that the vast majority of GNR samples are 
produced by lithographic patterning, characterized by 
rough edges at the atomic scaled"—. Already at low 
concentrations, such defects can destroy conductance 
quantizationii. Edge roughness may be largely sup- 
pressed in GNRs produced by unzipping single-wall car- 
bon nanotube.— However, the latter is not free of bulk 
defects. 

The way disorder affects electronic transport in GNRs 
strongly depends on its spatial range. For long-ranged 
disorder (LRD), corresponding to a ratio d/a^ ^ 1 be- 
tween the potential range d and the lattice parameter oq , 
backscattering is suppressed and the transmission is little 
affected. In contrast, short-ranged disorder (SRD) favors 
scattering processes with large momentum transfer such 
as backscattering. In this regime, quantum interference 
can cause wave function localization. In general, these 
simple arguments provide a qualitative explanation for 
the observed behavior of the conductance in current ex- 
periments. Edge roughness is essentially short ranged, 
while substrate impurity charges and ionically bonded 
adatoms are the typical sources of long-ranged disorder jii 

In view of the unavoidable disorder, the natural ques- 



tion that arises is whether one can indeed observe perfect 
transmission or conductance quantization in GNRs. This 
question has been theoretically investigated and partially 
answered by Wakabayashi et aZj ^^d^ . They have found 
that zigzag edge GNRs in the presence of long-range dis- 
order exhibits a quite robust perfectly conducting chan- 
nel (FCC). 

The dispersion relation for zigzag GNRs, Fig. (TJ helps 
one to understand the origin of FCC. Fig. [1] indicates 
two possible transversal momentum states for electronic 
energies such that, at low enough energies, only the first 
sub-band is allowed. The state close to the K point corre- 
sponds to a right propagating channel whereas the other, 
close to K' , is related to a left propagating channeli^. 
In this case, electrons propagate through the system in 
a well-defined direction leading to perfect transmission. 
Disorder can modify this scenario, provided it can cause 
a momentum transfer A/c \K — K'\ k, I/oq so that 
it mixes left and right propagating channels. In other 
words, this backscattering process requires a momentum 
transfer from states at the vicinity of the _ftr-point (rem- 
iniscent of bulk graphene) to the K' point, and vice- 
versa. The correspondence between short-ranged defects 
in rough-edged GNRs and intervalley scattering has been 
established by analyzing the scattering processes in both 
real and Fourier spaceii. For LRD such correspondence 
is more subtle and has been recently addressed by ana- 
lyzing the reflection probabilities of zigzag and armchair 
GNRs and their symmetry properties^. For disordered 
zigzag GNRs, electronic scattering should mix valleys, 
leading to the suppression of transmission depending on 
d. We address this issue by establishing a connection be- 
tween the FCC and the backscattering mechanism, ana- 
lyzing both the SRD and LRD regimes. 

The purpose of this work is to determine the funda- 
mental mechanisms that lead to the FCC by directly 
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Figure 1. Electronic bandstructure of a zigzag GNR with 
M = 10 horizontal chains along of width. K and K' are the 
two inequivalent points of the first Brillouin zone of graphene. 
The left (right) arrows indicate positive (negative) group ve- 
locity and forward (backwards) electronic propagation. 



analyzing the conductance dependence on the range of 
disorder as the scattering potential changes from short 
to long ranged. We have found that the PCC is not 
as robust as previous studies pointed out. In fact, we 
demonstrate that the emergence of the PCC crucially 
depends not only on the disorder range but also on the 
electronic energy. 

This paper is organized as follows. In Sec. Hi] we de- 
scribe the tight-binding model with on-site energy disor- 
der, the appearance of the PCC and the recursive Green's 
functions method. The numerical results are shown in 
Sec, mil where we discuss the effect of the transition from 
short- to long-ranged disorder on the conductance g (and 
also on the PCC). We also propose an analytical estima- 
tive for the degree of backscattering process to elucidate 
the physical origin of the PCC, which is compared to 
the conductance numerical results. Finally, Sec. IIVI is 
devoted to the conclusions. 



II. MODEL AND THEORY 

In this Section we present the model and the theory 
used to address the single-particle transport properties 
in GNRs. We first obtain analytical solutions for the 
band structure and wave functions of the tight-binding 
Hamiltonian that describes the electronic properties of 
pristine zigzag GNRs. Next, we introduce the local dis- 
order model used to numerically study the conductance 
in GNRs. We close the Section with a brief description 
of the numerical method employed to calculate the trans- 
port properties, namely, the recursive Green's functions 
method. 



A. Tight-binding model 

Close to half filling, the electrons in graphene are 
assumed to move by hopping through the Pz orbitals 
of the carbon atoms. Using the labels introduced in 
Fig. [21 the first-neighbor tight-binding Hamiltonian for 
graphene reads 

+ 4,m^«+l,m +H.C.^, (1) 

where the hopping parameter is t = 2.7 eV— and the sum 
is related to the sublattice A sites only. The operators 
al,m create and a„ „ (6„, annihilate an elec- 

tron at the site (n, m) of the sublattice A{B). The in- 
tegers n = 1, 2 and m = 1, . . . ,M label the atomic sites 
in the GNR unit cell according to the notation estab- 
lished in Fig. [51 M is related to the nanoribbon width by 
W = Maoy/3/2. The lattice parameter qq = 2.46 M re- 
lates to the carbon-carbon distance a through uq = ay/3. 




n = 1 2 ■•■ JV-1 JV N + 1 



Figure 2. (Color online) Zigzag graphene ribbon geometry, n 
and m label the atomic sites on the sublattices A (yellow dots) 
and B (blue dots). The atoms labeled n = 1, ■ • ■ , belong 
to the ribbon region and the atoms n = and n = -I- 1 
belong to the left and right contacts, respectively. 

The zigzag edge boundary conditions are^^ 

4,olO> = 0, (2) 
&Im+iIO>-0, (3) 

where |0) is the vacuum state. Notice that we require the 
states of each sublattice to vanish at opposite edges. 
The GNR eigenvalue problem reads 

H\a,k,,)=Eo,{k,)\a,k,), (4) 
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where is the longitudinal wavenumber and a the band 
(or channel) index. Using the on-site probability am- 
plitudes c^.n..mikx) = {0\an,Tn\a,k^) and cf.„„,(fc^) = 
(0| bn.m \a, k^) we write 



W, kx) — ^ (c^;n,m(^2:)^Q;n,m + l,m (^a: i^„i^ |0) 

(5) 

By inserting the eigenstates expansion (O into the 
Eq. ([1]), one obtains^S 



4„,m(fcx) = AA„(fc,)e™''°'===/2 sin[mz.,(fc,)] , 



(6) 
(7) 

where the minus and plus signs denote states with posi- 
tive and negative energies, respectively. The normaliza- 
tion factor 



AA„(fc,) = 1/ 



\ 



7V^sin2[ 



(8) 



is obtained by imposing (a, fc^-la, fc^:) = 1- The momen- 
tum function Va (kx) is introduced to satisfy the boundary 
condition Va{kx) is given by the multiple solutions 
of the transcendental equatio n^^i^° 



2 cos(aofc^/2) = — 



sm{Mva 



sin [{M + 



(9) 



These analytical solutions will be used in the calculation 
of the backscattering matrix elements. Finally, the zigzag 
GNR eigenergies read 



Ea{kx)lt = ± 



s\n[va{kx 



(10) 



B. Bulk disorder in GNR 

We calculate the conductance of disordered GNRs of 
length L = Nao/2 (Fig. [2]). To treat disorder, with em- 
ploy the Gaussian disorder model, defined as follows. We 
randomly choose Nij^p sites as the centers of Gaussian 
potentials with range d. A^imp is expressed in terms of 
the impurity concentration riimp = Mmp/-^tot, where 
the total number of atoms in the scattering region is 
^tot = NM. Hence, the disorder potential V at the 
position r reads 



l/(r) 



-R(Tii, mi)p/d^ 



(11) 



where R(ni, nii) is the center is the ith Gaussian disorder 
potential. The on-site lattice representation of V is 



V 



E 



— l,rn^n— l,m — 1 , m 



(12) 



where Vn^m corresponds to V{r) evaluated at R(n, m) 
corresponding to the position of the lattice site (n,m). 

The potential amplitude Ui is randomly chosen from a 
uniform distribution in the interval |J7i| < J7max, where 



/full space 



(13) 



The dimensionless parameter Uq/I defines the maximum 
disorder potential energy at each impurity site. 



C. Recursive Green's function technique 

The conductance is obtained by using the recursive 
Green's functions method^ii^^. This method provides 
a computationally efficient way to calculate the total 
Green's function of a GNR connected to pristine semi- 
infinite graphene leads at both ends. Using a decimation 
method we compute the surface Green's functions and 
the decay width functions of the left Fq and right Fjv+i 
leads. Next, we split the GNR "central" region (of width 
W and length L) in iV slices containing M transversal 
sites [see Figl2] and iteratively calculate the total re- 
tarded Green's function G\ jy that contains information 
about electron propagation from slice 1 to N . Finally, the 
dimensionless conductance g is obtained from the Caroli 
formula3^, namely, 5 = Tr \rQG\ j^TN+iG1j^] . The lin- 
ear electronic conductance is G = {2e^/h)g, where the 
factor 2 is due to the spin degeneracy. 



III. RESULTS 

In this Section we study the robustness of the PCC 
in disordered zigzag GNRs. This is done by numerically 
computing the dimensionless conductance g and inter- 
preting the results in terms of the analytical tools pre- 
sented in the previous Section. 

We compute the dimensionless conductance averaged 
over a large number of disorder realizations (g) (typi- 
cally 10'^) by means of the recursive Green's function 
method. The impurity potential strength is Uo/t ~ 1. 
As in Ref. HI, we consider GNRs with M = 10. For 
a zigzag GNR of this width, there is a single propagat- 
ing channel for E < E2 = 0.406i, where we denote the 
threshold energy to open the ath channel by Ea- 

We start by numerically investigating the behavior of 
the PCC in the crossover from SRD to LRD regimes as 
a function of the energy E. Figure [3] shows our results 
for (g) as a function of the potential range d/ao for two 
values of impurity concentration, namely, nimp = 0.10 
and riimp = 0.01. The simulations show that for suffi- 
ciently large values of d/ao, irrespective of E and n-m^p, 
the PCC always occurs, as (g) = 1 within the numeri- 
cal precision. Moreover, we find that the potential range 
dc, defined as the potential range above which the PCC 
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larger in the case of more diluted impurities, as expected. 
For both impurity concentrations, Fig. [3] shows conduc- 
tance plateaus for short scattering potential ranges, typ- 
ically d/ao < 0.3. These plateaus have a simple interpre- 
tation. Let us consider a single Gaussian disorder scatter- 
ing center, placed at a site i. For d smaller than roughly 
half the inter-atomic distance a/2 = ao/(2-\/3) = 0.29ao, 
the neighboring sites of i are hardly affected by the scat- 
tering center placed at i. Further reduction in the poten- 
tial range does not change the system Hamiltonian. 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 
E/t 



Figure 4. Disorder range above which the PCC appears, 
dc/ao, as a function of the electron energy E. Inset: Scaled 
average conductance {{g) — (7o)/(l — go) as a function of d/dc 
for different energies E/t. 



Figure 3. Average dimensionless conductance (g) as a func- 
tion of the potential range d for a ribbon of length L — 500ao 
and width W = 5%/3ao for the impurity concentrations (a) 
"■imp ~ 0.10 and (b) nimp = 0.01. 



appears, depends strongly and non-monotonically on en- 
ergy E: (i) dc decreases for increasing energies E start- 
ing from the charge neutrality point, and (ii) dc increases 
with increasing energy as E approaches £'2- This is in 
contrast to Ref. [l^ which suggests the emergence of 
the PCC for all energies E < E2 &t d/ao ~ 1-5 and 
nimp = 0.10. The same qualitative trend is found for 
both low (nimp = 0.01) and high (riimp = 0.10) impurity 
concentrations we analyze, as shown in Figs. [3^b) and 
^a), respectively. 

The overall values of the average conductance (g) are 



To investigate the role of the disorder potential range 
on the emergence of the PCC, in Fig.|3]we show dc/ao, de- 
fined above, as a function of energy E/t for nimp = 0.10. 
Figure [4] shows that dc/ao ~ 1, except for low energies 
and energies at the threshold of the n = 2 channel open- 
ing. The insert of Fig. 2] indicates that the behavior 
of the function {g{d/ao)) does not show, in general, a 
simple single-parameter scaling behavior. The numeri- 
cal results suggests that a single-parameter scaling for 
((5) ~ .9o)/(l ^ .90) ^ function of d/dc holds approxi- 
mately true only for 0.15 < E/t < 0.37, where go is the 
average conductance minimum for a fixed energy E/t. 

To understand the numerical results presented above 
we investigate the backscattering mechanisms induced by 
the disorder potential. These mechanisms can be quanti- 
fied by studying the backscattering matrix elements con- 
necting forward- to backward- moving states with a = 1, 
namely. 



{l,-k.\V\ 



(fc.)K 



-l;n — l,m 



i-k.)]'' 



(fcx)14- 



(14) 
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where the symbols are defined in Sec. |TI1 The expression 
for the backscattering matrix element (jl4p becomes very 
simple in the SRD regime: For a case of a single impurity 
placed at R„o,mo: it reads 

{l,-K{E)\V\l,+h:,{E)) = 

~ ['^rio.mo C-no ,mo{^x)Vno .mo , (15) 

where K = A or B. 

Let us examine the backscattering matrix elements in 
a number of representative situations. We first consider 
the single-impurity scattering case. Figure [5] shows the 
backscattering matrix, Eq. (fTSjl . as a function of the en- 
ergy for two different configurations: The Gaussian scat- 
tering potential is placed at the center or at the edge of 
the GNR. As for narrow zigzag GNRs and small values 
of E the electronic states are typically concentrated at 
the edges of the ribbon, we expect a distinct behavior in 
these two limiting cases. 

Figure [Sja) shows Vback(-B) = 
{l,—kx{E)\V \l,+kx{E)) as a function of d/ao for 
a single Gaussian potential placed at the center of the 
GNR. Due the SRD to LRD crossover, the backscat- 
tering matrix element decreases very fast as d/ag ^ 1. 
As expected from Eq. ([TC]) and the previous discussion, 
Vback(-E) hardly changes for d/ao < 0.25. For M = 10, 
the states \\,±kx) become localized at the GNR edges 
when E/t < 0.1—. In this situation, impurities located 
at the GNR center result in weak backscattering so 
that the plateaus quickly drop to 14ack(£') ~ 0. Figure 
\Elh) shows |Vback(^)| for the case of edge impurities. 
For E/t > 0.1, corresponding to ordinary states, the 
behavior is similar to Fig. [H^a). However, for edge 
states the situation changes dramatically: Here Vback(£') 
decreases surprisingly slowly with increasing d/ao > 1. 
This indicates that, in the presence of edge states, 
intervalley scattering is highly sensitive to the impurity 
position, even in the LRD case. 

We now consider the situations of low {rUmp = 0.01) 
and large (riimp = 0.10) impurity concentrations, the 
same values used in the conductance numerical calcu- 
lations (Fig. [3|). We compute the backscattering ma- 
trix element using Eq. (fT4|) and their average over dif- 
ferent configurations of disorder. The results are shown 
in Fig. [6] and arc contrasted with the average conduc- 
tance (g). Like (g), the average backscattering matrix 
elements, (|14ack(^^)P), also display a plateau-like be- 
havior for d/ao < 0.3. For \E\ > 0.15<, {\Vb^ck{E)\^) 
exhibits a fast decay with increasing d/aQ independent 
of the energy E. On the other hand, for \E\ < O.lSt 
the backscattering matrix elements decay nearly expo- 
nentially for 1.0 < d/aQ < 2.0, reaching minimum values 
depending on riimp. For energies near the Dirac point 
(E ~ 0), backscattering is maximal, persisting even in 
the LRD regime. In summary, in the SRD regime the 
averaged backscattering reaches its maximum value caus- 
ing a conductance suppression. In distinction, in general 
(|K>ack(£')p) is suppressed in the LRD favoring the ap- 
pearance of the PCC provided d/ao ^ 1- 
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Figure 5. Backscattering matrix element Vback calculated 
as a function of the scattering range d/ao for different values 
of the energy E for a single Gaussian impurity center placed 
(a) at the center of the GNR (no — 100, mo — 5) and (b) 
at its edge vicinity (no — 100, mo = 9). The GNR has the 
dimensions N = 200 and M = 10. For this value of GNR 
width, the edge state appears for E < 0.1^ .^^'^^ 



The small energy regime is an exception to this picture: 
The enhanced backscattering near the charge neutrality 
point can be understood in simple terms. Qualitatively, 
one expects that a scattering potential with a characteris- 
tic length scale d can only effectively backscatter electron 
states with initial momentum k to final momentum —k 
provided l/d w |— 2fc-|-G|; where G is a reciprocal lattice 
vector. Nonzero G vectors describe Umklapp scatter- 
ing processes. As one approaches the charge neutrality 
point, scattering between k « n/a to — fc « — 7r/a medi- 
ated by a G = 27r/a Umklapp becomes dominant, so that 
|— 2fc -|- G| approaches zero and even a very smooth po- 
tential (large d) is able to scatter the electron states and 
destroys PCC. That is the reason for the sharp increase 
of dc near E = 0, shown in Fig. |4l 

We discuss now the robustness of the PCC for ener- 
gies E approaching E2. As discussed, for such energies 
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Figure 6. Average backscattering matrix element (|VbackP) 
as a function of the disorder range d/ao for a number of en- 
ergies E. The average is taken over 10^ disorder configura- 
tions for an impurity concentration of (a) riimp = 0.01 and 



(b) 



= 0.10. The GNR has the dimensions L = 500ao 



and W — 10(\/3/2)ao. This ribbon presents edge state for 
E < O.lt. 



Vback ~ as long as d/ao > 1. To explain deviations from 
the PCC, first-order perturbation theory is not sufficient 
since LRD cannot account for the large momentum trans- 
fer necessary for backscattering, as already discussed. 
LRD can only account for backscattering processes in 
the vicinities of the energy E2 , where left and right prop- 
agating modes (Fig. [1]) are not far away in momentum 
space. For a sufficiently smooth long-range disorder, one 
can account for disorder effects by introducing a local 
chemical potential as /.t[y(r)]. In this scenario, LRD can 



suppress the PCC at energies E2 ~ SE < E < E2 + SE, 
where SE is the typical disorder potential fluctuation. 

We have also carried out numerical simulations for 
higher energies (up to E/t = 0.9), corresponding to the 
multi-channel case. However, in this case we have not 
found any qualitative difference between our results and 
the ones reported in Ref. [TgI . More precisely, we have 
found deviations from the PCC for incident energies at 
the vicinity of energies corresponding to the crossover be- 
tween g = 2n — 1 and 2n + 1 for clean ribbons (n being 
the number of channels) . This behavior corroborates the 
results of Ref. 



IV. CONCLUSIONS 

In summary, the electronic conductance of disordered 
GNRs in general, and the emergence of the PCC in par- 
ticular, present a richer behavior than expected from pre- 
vious studies. Specifically, the critical impurity poten- 
tial dc for the PCC emergence shows a strong and clear 
dependence on the electronic energy. This result, ob- 
tained numerically from recursive Green's functions cal- 
culations, is explained in simple physical grounds by cal- 
culations of energy and potential range dependences of 
electron inter valley scattering probabilities, which follows 
the opposite qualitative trends of the conductance. This 
occurs for sufficiently low energies such that the Fermi 
level crosses just one band (single channel conductance). 
This behavior confirms and justifies the simple picture of 
an electron undergoing intervalley scattering only if the 
impurity potential has substantial Fourier components at 
large momenta, thus being able to provide a momentum 
transfer Ak 1/ao to the electron. In other words, for 
a single channel, backscattering occurs only if interval- 
ley scattering does. This picture not only explains the 
conductance behavior for low energies, but it also pro- 
vides us with guidelines to explain the more complicated 
cases in which the Fermi level crosses many channels. In 
that case, as shown in Fig. 1, backscattering (i.e., group 
velocity reversal) can occur even for small momentum 
transfers (intravalley scattering), therefore explaining the 
disappearance of the PCC as the Fermi level approaches 
the threshold for opening the second transmission chan- 
nel. 
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